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ABSTRACT 

■ Two disks of young stars have recently been discovered in the Galactic Center. 

The disks are rotating in the gravitational field of the central black hole at radii r ~ 
£\j . 0.1 — 0.3 pc and thus open a new opportunity to measure the central mass. We find 

that the observed motion of stars in the clockwise disk implies M = (4.3 ±0.5) x 1O 6 M0 
for the fiducial distance to the Galactic Center Rq = 8 kpc and derive the scaling of 
t""". \ M with R . As a tool for our estimate we use orbital roulette, a recently developed 

_ ' method. The method reconstructs the three-dimensional orbits of the disk stars and 

J> ■ checks the randomness of their orbital phases. We also estimate the three-dimensional 

m 

(N 



positions and orbital eccentricities of the clockwise-disk stars. 
Subject headings: Galaxy: center — galaxies: general — stellar dynamics 

1. INTRODUCTION 



o 

\o 
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A supermassive black hole is believed to reside in the dynamical center of the Galaxy (Eckart 
& Genzel 1996; Genzel et al. 1997; Ghez et al. 1998; Schodel et al. 2002; Genzel et al. 2000; 
Ghez et al. 2000; Schodel et al. 2003; Ghez et al. 2003). It is associated with the compact radio 
source Sgr A* and dominates the gravitational potential at distances r < 1 pc from the center. The 
^ . central mass has been estimated with a few different methods. In particular, the statistical analysis 

of sky-projected positions and velocities of stars in the central 0.5 pc gave estimates that range 
above 3 x 1O 6 M (Genzel et al. 2000; Schodel et al. 2003). The systematic error of these estimates, 
however, is difficult to quantify as they depend on the assumed three-dimensional statistical model 
of the cluster. 

Independent estimates of the black hole mass were derived from observations of a few stars 
(S-stars) that move on very tight orbits around Sgr A* , at radii r ~ 0.01 pc, and have orbital 



1 Physics Department and Columbia Astrophysics Laboratory, Columbia University, 538 West 120th Street New 
York, NY 10027 

2 Astro-Space Center of Lebedev Physical Institute, Profsojuznaja 84/32, Moscow 117810, Russia 
3 Leiden Observatory, PO Box 9513, NL-2300 RA, Leiden, The Netherlands 
4 Max-Planck Institut fur extraterrestrische Physik (MPE), Garching, Germany 
5 Department of Physics, University of California, Berkeley, CA 94720 



- 2 - 



periods as short as a few decades. Fractions of S-star orbits have been mapped out over the past 
decade (Schodel et al. 2002; Schodel et al. 2003; Ghez et al. 2003; Ghez et al. 2005). By checking 
which value of the central mass provides an acceptable fit to the orbital data, Eisenhauer et al. 
(2005) found M = (4.1 ± 0.4) x 1O 6 M , and Ghez et al. (2005) found M = (3.7 ± 0.2) x 1O 6 M 
for the fiducial distance to the Galactic Center R = 8 kpc. Same star orbits were used to derive 
Ro = 7.62 ± 0.32 kpc (Eisenhauer et al. 2005), in agreement with previous methods (e.g. Reid et 
al. 1993). 

In this paper, we report a new independent estimate of the central mass that uses the young 
massive stars at r ~ 0.1 — 0.3 pc as test particles. A remarkable property of these stars has been 
discovered recently: they form a disk-like population (Levin & Beloborodov 2003; Genzel et al. 
2003). Two stellar disks apparently coexist at these radii. One of them rotates clockwise on the 
sky and the other one — counter-clockwise. New observations of Paumard et al. (2006) have 
increased the number of known disk members by a factor of three and show clearly the two-disk 
structure of the young population. The clockwise disk is especially well pronounced. It has about 
30 members and a small thickness H/r « 0.1. 

In this paper, we investigate the orbital motions of stars in the clockwise disk and infer the 
central gravitating mass. As a tool of our mass estimate we use orbital roulette, a recently developed 
statistical method (Beloborodov & Levin 2004, hereafter BL04). This method is generally more 
precise than the estimators based on the virial theorem and provides a statistically well-defined 
way to quantify the uncertainty of the mass estimate. 

2. CLOCKWISE DISK 

The disk-like structure of the young stellar population is not evident when it is observed in 
projection: the line-of-sight coordinaties z of the stars are unknown. This structure, however, is 
evident in the distribution of the measured 3D velocities. In particular, the velocity vectors of stars 
that move clockwise on the sky cluster near a particular plane. This may be only if the orbits of 
stars are near this plane and form a disk-like population. 

The clockwise disk has a small thickness at angular distances 0.7" < d < 5" from Sgr A* . 
We use the stars in this region to define the midplane of the disk. It is defined as the midplane 
of the observed velocity vectors, and its normal vector is 1 no = (—0.135,-0.857,0.497). The 
deviations of individual stellar velocities from the midplane are shown in Figure 1. Hereafter in 
this paper we use the sample of stars between 0.7" and 5" with firmly detected clockwise rotation 
l z = xv y — yv x > 2Al z , where Al z is the measurement error in l z . 



1 (x,y,z) are standard coordinates: they are chosen so that Sgr A* is at the origin, the x and y axes are in the 
plane of the sky and directed along increasing right ascension and declination, respectively, and the z axis is directed 
along the line of sight away from us. 
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Fig. 1. — Quantity \i = no • v/v for all young stars observed by Paumard et al. (2006) that have 
apparent clockwise rotation on the sky l z > 0, shown versus projected distance from Sgr A* . For 
randomly oriented orbits, fi would be uniformly distributed between —1 and 1. The actual data 
show a strong clustering around [i = (midplane). Data shown in black represent the stars with 
firmly detected clockwise rotation l z > 2AZ 2 , and cyan points — the stars with < l z < 2Al z , 
where Al z is the error in l z . The sample used in this paper is the 28 stars shown in black in the 
region 0.7" < r± < 5". 



-4 - 



2 - 







"i 1 1 r 



"i 1 1 r 



"i 1 1 r 







J I I L 



J I I L 



I I I L 



n 1 1 r 



N = 28 

< f\> = 7° 



i — i . 



10 20 30 

deviation [degree] 



40 



Fig. 2. — Histogram of the angles ipi betwen observed best-fit velocities Vj and the disk midplane 
with normal no- Angle ip is related to \x shown in Figure 1 by sin-*/? = [i. The mean value of |^>| 
is (|V>|) = 7.1° and (tp 2 ) 1 / 2 = 9.3°. Dashed curve shows the distribution of \ip\ expected from our 
model disk population. The dispersion of orbital inclinations 6 in the model is 10°. 
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2.1. Thickness of the Disk 

The sample has 28 stars, which is sufficiently large to study the distribution of velocity devia- 
tions from the disk midplane (Fig. 2). The deviation angles ipi are defined by ipi = arcsin(no • Vi/v «). 
The observed distribution of tp m ay be fitted by a disk population model. We will use a simple 
Gaussian model of orbital inclinations 9 around the midplane. The probability distribution of 9 is 
then given by P{9) = (2/vr) 1 / 2 (A^)- 1 exp[-(^/A0) 2 /2]. 

The Gaussian model of the disk has one parameter — the dispersion A9 — and predicts a 
certain probability distribution for -ip. We calculate this distribution using Monte-Carlo technique, 
including a simulation of observational errors. We randomly draw N = 28 orbits from the model 
population and choose a randomly directed vector of length Vi in each orbital plane i = 1,...,N. 
This gives 28 simulated velocities Vj. Then observational errors Avj are added to each Vj (the 
real vi and Avj from Paumard et al. 2006 are used). Finally, deviations ^ of velocities Vj from 
the disk midplane are calculated. This simulation is repeated 10 5 times, which gives the accurate 
distribution of ip predicted by the model as well as the average values (1-01) and (tp 2 ) 1 / 2 as functions 
of A9. 

The observed (|^|) « 7° and (ip 2 ) 1 / 2 ss 9° are reproduced by the Monte-Carlo model with 
A9 « 10°. The model is consistent with the observed distribution of tp for A9 = 10 ± 3°. Hereafter 
we will use the best-fit model with A9 = 10°. 

The geometrical thickness of the disk may be described by the elevation of stars above/below 
the midplane, 

— = I sin <p sin 01, (1) 

r 

where (j) is the orbital phase of the star, and 9 is the inclination of its orbit. Assuming circular 
orbits for simplicity, one finds after averaging over the orbital phase, 

7 = (-J < sin *>- ( 2 ) 

This gives for the Gaussian model (approximating sin(#) « (sin9) for small A9 = 0.17), 

H ( 2 \ 3/2 

— = - sin A9 = 0.09 ± 0.03. (3) 
r \ir J 



2.2. Circular Motion? 

We now check whether the data is consistent with circular motion of the disk stars. The 
assumption of circular motion implies that radius-vector rj of the i-th star satisfies 

rf r -v 4 = 0, i = l,...,N. (4) 
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This equation, together with the known position (xj,yj) on the sky, determines Z{ — the line-of- 
sight coordinate of the star. The easiest way to test the circular-motion hypothesis is to look at 
the orientations of orbital planes implied by this hypothesis. The inferred rf r allow us to calculate 
the angular momentum vector for each star 

lf = rfxv,, i = l,...,N, (5) 

and then find its deviation 6f r from the disk axis no- The dispersion of 6f T would be expected 
to be consistent with disperion A9 « 10° of orbital planes inferred directly from Vj. The actual 
distribution of 9f T is shown in Figure 3. One can see that the circular hypothesis implies orbital 
inclinations that are inconsistent with the found thickness of the disk. A significant fraction of the 
disk stars must have non-negligible eccentricities. 

A circular orbit would provide a straightforward estimate of the central mass M. Indeed, 
circular motion implies a simple relation between M, v, and r: GM = rv 2 . An estimate of M 
could be done by applying this relation to stars with small 8 ai , which are possibly on circular 
orbits. However, it is difficult to quantify the accuracy of such analysis. A more accurate estimate 
of M is made in § 4, without making any prior assumptions on eccentricities. 



2.3. Three-Dimensional Positions 

Measurements of stellar positions relative to Sgr A* are limited to 2 dimensions — the plane 
of the sky. However, the stars that belong to the clockwise disk have an extra constraint on their 
3D positions r^: they must be close to the midplane of the stellar disk. This constraint may be 
written as 

(n ± An) • r, = 0, i = l,..,N, (6) 

where no is the disk axis and An represents the dispersion of orbital planes around the disk 
midplane. 

The line-of-sight coordinate of the i-ih star is then given by 

Zi = ■ — , nj = n ± An 

The orbital plane nj may differ from the best-fit disk plane by ~ 10°, which induces an error in 
the inferred z-coordinate. To estimate the error we rock nj with a Gaussian distribution around 
no with a mean deviation of 10°. 

The resulting Zj and their uncertainties are found using the following Monte-Carlo simulation. 
We draw randomly nj and v x i, v y i from their measured Gassian distributions. Then we determine 
v Z j from the condition Vj ■ nj = 0. The obtained realization of nj and Vj is given a weight that 
corresponds to the level of consistency of v z> i with the measured Gaussian v Z j = v z a ± Aw Zj j. (In 
the simulation, we get this weight by rejecting the realization if its v z ^ is farther in the Gaussian 
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Fig. 3. — Histogram of orbital inclinations 9 that would be inferred from the circular hypothesis. 
Dashed curve shows the 10° Gaussian distribution of 9 that fits the 3D velocity data (see Fig. 2). 
The circular hypothesis implies large orbital inclinations 9 > 30° for at least 1/4 of stars, which is 
inconsistent with the observed narrow distribution of velocity angles ipi- 
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tail than a randomly chosen v' zi from the Gaussian distribution.) In this way, a realization rij, Vj 
is given the weight that takes into account both the thickness of the disk and the error of the 3D 
velocity. We find the mean Zi and its dispersion Azi using about 10 6 realizations. 

The deduced 3D positions together with the measured 3D velocities give 6D information for 
each star. Then the velocities v r and may be deduced, where r and cf> are polar coordinates in 
the orbital plane of the star. The measurement errors of Vj and the uncertainty in the individual 
orbital planes imply errors Av r and Av^. We estimate the uncertainties assuming the Gaussian 
dispersion in rij with AO = 10° and Gaussian velocity errors Avj. The results are summarized in 
Table 1 and illustrated in Figure 4. 



3. MASS ESTIMATION USING ORBITAL ROULETTE 

The orbital periods of stars at r ^ 0.1 pc are thousands of years, and it is presently impossible 
to trace any significant fraction of the orbits. Therefore, a statistical method must be applied using 
the observed instantaneous positions rj and velocities Vj of the stars. Such a method must take 
into account the possible errors in rj and Vj. 

A traditional mass estimate would be based on the virial theorem, 

The virial relation gives M = 3.8 x 10 6 Mq. One faces difficulties, however, when trying to define 
the error of this estimate (see BL04 for discussion). 

A major assumption in the derivation of the virial relation (8) is that the stars have random 
orbital phases. BL04 have shown that the random-phase ("roulette") principle can be applied 
directly to data, without invoking the virial theorem, and this gives a more accurate and statistically 
complete estimator. It judges any trial gravitational potential 3>(r) by reconstructing the orbits in 
this potential and checking the inferred orbital time phases. The true potential must give a "fair 
roulette" — a random (flat) distribution of the phases. 

The 3D positions of the disk stars have uncertainties because their precise orbital planes 
are unknown. Besides, the 3D velocities are measured with errors. We take into account the 
uncertainties using Monte-Carlo technique: many random realizations of the data set (rj,Vj) are 
generated and studied. The possible (rj, Vj) are generated assuming the 10° Gaussian distribution 
of the orbital planes and Gaussian errors in Vj with dispersions Avj given in Paumard et al. (2006). 

Each realization of the data set (rj,Vj) is analyzed by the roulette estimator as described in 
BL04. The analysis aims to accept or rule out a trial potential 5>(r) = — GM tv \ & \/r at a chosen 
confidence level. Given M tria i and a data set (rj, Vj) one can reconstruct all N = 28 orbits and find 
the time orbital phases gi of the stars. The time phase g is defined so that g = at the pericenter 
of the orbit and g = 1 at the apocenter. To judge M tr i a i we check the consistency of with the 
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Fig. 4. — Positions and velocities of the 28 stars in the disk plane (p, q) . Sgr A* is located at 
the origin (black circle). The p-axis is chosen in the intersection of the disk plane with the (x,z) 
plane and the g-axis is perpendicular to the p-axis. The axis directions (unit vectors) e p and e q 
are chosen so that e p ■ e x > and e q ■ e y > 0. The velocity of each star is shown by 5 vectors: the 
best-fit vector Vj, Vj ± Av p e p , and Vj ± Av q e q . We show the stars in a single plane, as if they had 
a common orbital plane. Their real orbits may differ from the disk midplane by ~ 10°, and the 
corresponding uncertainties in positions are indicated in the figure. The distance to the Galactic 
Center Rq = 8 kpc is assumed. 
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expected flat distribution between and 1 (and do that for many realizations of the data sets to 
take into account the uncertainties in the data). 

The roulette algorithm is described in detail in BL04. Here we use its simplest version based 
on the analysis of the mean phase g = N^ 1 J2 9i- The true M must satisfy the condition 

g(M) = 0.5 ± (12iV)~ 1/2 . (9) 

If M tr iai is too small then g{M tT - ia [) — ► 0, i.e. all stars will be found near the pericenter, while too 
large M tria i will give g(M tria i) — > 1, i.e. all stars will sit near apocenter. One can accept only such 
Atrial that give g consistent with the fair roulette (eq. 9). 

In our Monte-Carlo simulation, <?(M tr j a i) that is found for a realization of the data set (rj, Vj) 
is compared with g' found for a realization of a fair roulette g[. The statistics of this comparison is 
accumulated for 10 6 random realizations. Thus, the probability P_ of g(M tria i) < g' is accurately 
found. 

The probability P_(M tr i a i) is monotonically increasing with M tr i a i since larger masses imply 
larger orbital phases (closer to the apocenter). Trial masses that give small P_ are ruled out with 
significance P_ (BL04). If P_ > 1/2, we look at P + = 1 — P_. Trial masses that give small 
P + < 1/2 are ruled out with significance P+. 

The found functions P±(M t n a \) are shown in Figure 5. The confidence intervals for M are 
between the two curves. For example, the 80% confidence interval (P_ = P+ = 0.1) is 3.85 x 
1O 6 M < M < 4.93 x 1O 6 M0. The confidence intervals approximately correspond to a Gaussian 
distribution with standard deviation (0.4 — 0.5) x 1O 6 M0, and we can summarize the result as 
M = (4.3toj) x 1O 6 M . 

We have also checked the distribution of gi for the found M. It is consistent with the expected 
flat distribution; this is verified using the Anderson-Darling measure (see BL04). We conclude 
that the observed stellar motions are consistent with Kepler potential = —GM/r with M = 
(4.3183) x 1O 6 M . 

There are two sources of systematic error of our result. First, we had to chose a disk-population 
model in the calculations. We specified the midplane (normal no) and the Gaussian distribution of 
orbital inclinations around the midplane with dispersion A9. The parameters no and AO have been 
derived by fitting the 3D velocity data and have finite errors. The obtained M is most sensitive 
to n Z Q which determines the disk inclination to the line of sight. Our best-fit disk has n Z Q ~ 0.5; 
larger n z would imply smaller sizes of the orbits and hence smaller M. n Z Q as large as 0.6 may 
still be consistent with the velocity data (Paumard et al. 2006), and it would change our estimate 
to M = (4.0lg;|) x 1O 6 M . The parameter A6 has the error of 3° - 4°. The choice of AO = 13° 
instead of 10° would broaden the confidence intervals by a factor ~ 1.2. 

The second possible source of error is the distance to the Galactic Center Rq which we assumed 
to be 8 kpc. The mass estimate depends on Rq because the proper velocities vj_ = (v x ,v y ) and 
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positions (x, y) scale linearly with Rq (v z is not affected as its measurement is based on the Doppler 
effect). The midplane of the disk varies little and may be assumed constant to a first approximation. 
Then the deduced radial positions r of the stars scale linearly with Rq. The dependence of M on 
Rq may be derived using a toy problem with N 3> 1 stars moving in the disk midplane on a circular 
orbit of radius r true . The true and deduced masses are then given by 

GM tmc = r truc (wi )true + vl) , 



GM = r (vl + vl) = r truc (-^-\ 

\ -Kotruc / 



y ±,true 



Rq 



-Rotrue / 



After averaging over N stars and taking into account that (v± )trnc /v tIUC ) = (l+n 2 z )/2 and (v z /vt IVS ) = 
(l-n 2 z )/2, we find 



M(i? truc) = M(8 kpc) 



2< 1 + n '\ 8 kpc 



Ro 



true 



-Rotruc 

8 kpc 



(10) 



We have empirically checked this dependence by repeating the roulette simulation for various Rq 
(as small as 7 kpc). The obtained P±(M) are well described by a simple shift of M according to 
equation (10). In particular, Rq = 7.6 kpc (Eisenhauer et al. 2005) would give M = (3.8!q 3 ) x 
1O 6 M . 

Table 1 shows the eccentricities of reconstructed star orbits for Rq = 8 kpc and the best-fit 
central mass M = 4.3 x 1O 6 M . The errors Ae^ [i = 1, ...,28) are caused by uncertainties Avj and 
Arij. The errors are found using the Monte-Carlo simulation that is similar to the estimation of 
Azi (see § 2.3). This simulation gives the probability distribution /(e«) for each e^. A non-zero 
eccentricity is detected if > 2Ae« and accurately measured if > 3Aej, where is the mean 
of /(ej) and Ae^ is its dispersion. If distribution /(e^) is broad (large errors), its mean represents 
the errors rather than the true e^. Therefore, e« is an accurate estimate of eccentricity only if 
the distribution /(e^) is sufficiently narrow, Ae« ^ 0.3ej. Table 1 shows e,- L ± Ae« for stars with 
Aej < 0.5ej. For stars with Ae^ > 0.5ej we give + Ae^ as an upper limit. The obtained values 
are consistent with the results of Paumard et al. (2006) where eccentricities are estimated with 
another method. 



4. DISCUSSION 

The reconstructed orbits of stars in the clockwise disk depend on the assumed central mass 
M. In particular, the star orbital phases depend on M. We have used this dependence to estimate 
M: the true M must give time phases consistent with a uniform distribution. This is the first 
practical application of orbital roulette, a new method of estimating the gravitational potential 
from instantaneous motions of test bodies (BL04). 
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Fig. 5. — Confidence intervals for Sgr A* mass M found with the roulette method. For example the 
80% confidence interval (P+ = P_ = 0.1) is shown by the horizontal solid line (P+ = P_ = 0.1). 
The distance to the Galactic Center Rq = 8 kpc is assumed. The scaling of the estimated mass 
with Rq is given by eq. (10). 
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The obtained estimate M = (4.3^04) x 10 6 M Q is consistent with the recent estimates based 
on the maps of orbits of S stars in the very vicinity of the center (r < 0.01 pc). The independent 
estimate presented here uses different stars at r = 0.1 — 0.3 pc and is complementary to the small- 
scale estimates. It measures the total mass within r ~ 0.1 pc which may include stars and dark 
matter. 

The approximate equality of the masses within the central 0.1 pc and 0.01 pc implies an upper 
bound on the mass between 0.01 and 0.1 pc: AM & 0.5 x 1O 6 M . This constrains a possible cusp 
of the stellar population around the cenral black hole as well as the dark-matter content of the 
central region. Further observations may provide a better constraint or resolve AM. 

The data shows an inner cutoff of the disk population at a radius r ~ 0.05 pc (see Fig. 1 and 
Paumard et al. 2006). The observed peak of the population at r ~ 0.1 pc supports the hypothesis 
that the young stars were born in a gravitationally unstable accretion disk (LB03). 

The found geometrical thickness of the clockwise stellar population H/r = 0.09 ± 0.03 (eq. 3) 
does not imply that it was born with this H/r. The parent accretion disk was likely thinner, and 
its stellar remnant have thickened with time. The age of the stellar population, t age = (6 ± 2) 
million years, should be compared with the timescale of diffusion of orbital parameters. The latter 
may be written as 

to ~°-Wp*m(iV2)' 

where m* is a characteristic mass of the objects that dominate fluctuations of the gravitational 
potential, p* is their mass density averaged over the considered region r ~ 0.3 pc, and N* ~ 
(4ir/3)r s p :¥ /m > t is the total number of these objects. In a uniform stellar core of total mass ~ 
3 x 10 5 M Q , the fluctuations would be dominated by low-mass stars, whose total number is iV* ~ 10 6 . 
Substituting the observed characteristic velocity, v ~ 300 km/s, one would find tdiff ~ 10 10 yr3> 

tage- 

However, the Galactic Center is not uniform. The two stellar disks and stellar clusters are ob- 
served, which significantly perturb the gravitational potential. These perturbations should strongly 
reduce For example, if 1/30 of the core mass resides in iV* ~ 10 objects with m* ~ 1O 3 M then 

tdi/f ~ 10 8 yr. The corresponding acquired thickness of the stellar disk is H/r ~ t age /tdiff ~ 0.1 in 
agreement with the observed value. A similar estimate for tdiff may be derived from the interaction 
between the two young disks (Nayakshin et al. 2005). 

The eccentricities of disk stars must also diffuse with time. If the stars were born on circular 
orbits then e ~ 0.1 are expected to develop as the disk thickens to H/r ~ 0.1. Table 1 shows that 
a significant fraction of the young stars in the clockwise disk have relatively large eccentricities. 
The data are at best marginally consistent with the hypothesis that the stars used to have circular 
orbits 6 million years ago, and we leave a detailed analysis to a future work. Note that the stars 
may have been born in a gaseous disk that was not completely circularized: the circularization time 
was longer than the timescale for star formation. Then the initial eccentricities e of the stars are 
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non-zero. Besides, the interaction of stars with the disk might have influenced e before the disk 
was gone. 

An alternative scenario assumes that the stars were originally bound to an inspiralling intermediate- 
mass black hole (IMBH; Hansen &; Milosavljevic 2003), so that the observed disk plane is associated 
with the orbital plane of the IMBH. In this scenario, the inspiralling young cluster was stripped by 
the tidal field of Sgr A* and left a disk of stars with a small dispersion of orbital planes. Calcula- 
tions of Levin, Wu, & Thommes (2005) show that the stripped stars can generally have eccentric 
orbits. However, the inspiralling cluster scenario appears to be in conflict with the observed stellar 
distribution in the clockwise disk (see also Paumard et al. 2006 for discussion). 

AMB was supported by Alfred P. Sloan Fellowship. 
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Table 1: Positions and velocities of the clockwise disk stars in their orbital planes 



Name 


i> 


X 


y 


z 


r 


4> 






e 


E15: [GKM98] Sl-3 


0.5 ±2.9 


0.52 


1.03 


1.88 ± 0.40 


2.22 ± 0.34 


80 ±3 


-12 ±42 


534 ± 22 


< 0.33 


E16: [GEQ97] W5 


13.7 ±4.7 


-1.13 


0.31 


0.27 ±0.28 


1.23 ±0.09 


164 ± 16 


51 ± 187 


608 ± 63 


< 0.47 


E17 




-6.9 ±3.8 


-0.05 


-1.21 


-2.04 ±0.57 


2.40 ±0.52 


271 ± 2 


-33 ± 34 


429 ± 29 


0.28 ±0.13 


E18: 


[GE097] Wll 


-11.4 ±3.4 


-0.79 


-1.04 


-2.35 ± 0.67 


2.71 ±0.59 


251 ± 7 


274 ± 68 


408 ± 45 


0.62 ±0.14 


E20: 


GCIRS 16C 


-18.5 ±4.9 


1.36 


0.58 


1.47 ±0.56 


2.11 ±0.42 


59 ± 9 


-70 ± 83 


454 ± 48 


0.34 ±0.14 


E21: 


[GE097] W13 


-1.4 ±4.0 


-1.02 


-1.20 


-2.29 ± 0.54 


2.81 ±0.50 


252 ±4 


-99 ± 38 


389 ± 29 


0.35 ±0.11 


E22: 


[GE097] W10 


0.2 ±3.7 


-1.63 


-0.37 


-1.08 ± 0.40 


2.02 ±0.22 


217 ± 11 


148 ± 106 


505 ± 45 


0.35 ±0.15 


E23: 


GCIRS 16SW 


6.4 ±5.8 


1.26 


-1.18 


-1.75 ± 0.61 


2.50 ±0.48 


308 ± 12 


-131 ±85 


398 ± 49 


0.42 ±0.16 


E24: 


[GE097] W7 


-0.2 ± 5.0 


-2.00 


0.17 


-0.25 ± 0.46 


2.08 ±0.13 


188 ± 14 


-36 ± 99 


390 ± 43 


0.46 ±0.12 


E25: 


[GE097] W14 


-12.3 ± 5.1 


-1.97 


-0.60 


-1.38 ± 0.79 


2.58 ±0.57 


216 ± 15 


-86 ± 94 


329 ± 49 


0.54 ±0.15 


E27: 


GCIRS 16CC 


-9.3 ±6.5 


2.41 


0.65 


1.62 ± 0.65 


3.03 ±0.39 


40 ± 12 


95 ± 79 


300 ± 43 


0.54 ±0.15 


E28: 


GCIRS 16SSE2 


0.0 ±3.4 


1.74 


-1.79 


-2.55 ± 0.65 


3.61 ±0.48 


305 ±9 


-119 ±67 


404 ± 31 


0.37 ±0.14 


E29 




-14.2 ±6.6 


1.19 


2.20 


3.85 ± 0.68 


5.25 ±1.36 


85 ±4 


-97 ±49 


252 ± 30 


0.51 ±0.13 


E30: 


GCIRS 16SSE1 


-4.7 ±3.9 


1.91 


-1.63 


-2.22 ± 0.69 


3.41 ±0.52 


310 ± 11 


-27 ±69 


370 ± 28 


0.27 ±0.13 


E32: 


MPE+1.6-6.8 


6.4 ±6.7 


2.22 


-1.38 


-1.82 ± 0.76 


3.25 ±0.47 


321 ± 11 


-121 ± 101 


387 ± 56 


0.42 ±0.18 


E34:MPE+1. 0-7.4 


1.7 ±7.1 


1.52 


-2.26 


-3.24 ±0.73 


4.40 ±0.81 


297 ±7 


27 ±47 


313 ±43 


0.33 ±0.16 


E35: 


GCIRS 29NE1 


-2.9 ± 7.6 


-1.19 


2.47 


3.59 ± 0.72 


4.68 ±0.75 


111 ±6 


29 ±68 


372 ± 48 


< 0.59 


E36 




-1.7 ±4.4 


0.54 


2.75 


4.23 ± 0.68 


5.54 ±0.92 


89 ± 2 


45 ±23 


326 ± 29 


< 0.38 


E38 




-1.4 ±6.5 


0.23 


3.31 


4.16 ± 0.64 


6.35 ±0.98 


93 ± 3 


60 ± 25 


336 ± 40 


< 0.68 


E40: 


GCIRS 16SE2 


-3.7 ±9.4 


3.53 


-1.43 


-1.44 ±0.85 


4.15 ±0.33 


335 ± 11 


-82 ±95 


365 ± 53 


0.40 ±0.19 


E41: 


GCIRS 33E 


16.0 ±9.4 


0.78 


-3.74 


-4.29 ± 0.60 


7.85 ±1.65 


283 ±6 


-143 ±32 


214 ±47 


0.61 ±0.14 


E43 




-19.8 ± 7.2 


-1.92 


-3.35 


-4.07 ±0.71 


7.45 ±1.82 


253 ± 7 


16 ±55 


238 ± 32 


0.38 ±0.17 


E44 




-13.3 ±5.9 


1.74 


3.54 


4.36 ± 0.51 


8.14 ±1.58 


86 ±4 


-124 ±40 


252 ± 30 


0.53 ±0.16 


E50: 


GCIRS 16SE3 


-5.5 ±4.1 


4.02 


-1.39 


-1.12 ± 0.91 


4.49 ±0.26 


340 ± 10 


-121 ±83 


308 ± 37 


0.41 ±0.18 


E54 




12.6 ±5.0 


-4.40 


2.16 


2.55 ± 1.05 


5.75 ±0.76 


144 ± 13 


51 ± 72 


282 ± 31 


0.33 ±0.16 


E56: 


GCIRS 34W 


1.5 ±4.7 


-4.86 


1.91 


1.86 ± 1.06 


5.66 ±0.43 


157 ± 12 


-83 ± 79 


323 ± 34 


0.38 ±0.16 


E57 




3.0 ±6.9 


5.30 


0.30 


1.78 ± 1.15 


5.79 ±0.62 


24 ±11 


-33 ± 65 


240 ± 35 


0.45 ±0.16 


E61: 


GCIRS 34NW 


10.3 ±5.5 


-4.48 


3.42 


3.70 ± 0.97 


7.33 ±1.04 


132 ± 10 


-4 ± 63 


283 ± 29 


< 0.46 



Positions are offsets from Sgr A* . 

x, y, z, r are in units of 10 17 cm, v r and — in km/s, ip and <j) — in degrees. 

In order to express distances in arcseconds one should divide x, y, z, r in the table by 1.20. 



